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We describe the current status of the Ca optical frequency standards with laser-cooled neutral atoms realized in two 
different laboratories for the purpose of developing a possible future optical atomic clock. Frequency measurements 
performed at the Physikalisch-Technische Bundesanstalt (PTB) and the National Institute of Standards and 
Technology (NIST) make the frequency of the clock transition of ^°Ca one of the best known optical frequencies 
(relative uncertainty 1.2 ■ 10~^*) and the measurements of this frequency in both laboratories agree to well within 
their respective uncertainties. Prospects for improvement by orders of magnitude in the relative uncertainty of 
the standard look feasible. 



Nous decrivons le present etat de I'art des etalons de frequence optique avec des atomes neutres de Ca refroidis 
par laser, realises dans deux differents laboratoires, dans le but de developper eventuellement une future horloge 
atomique dans le domaine optique. Les mesures de frequences realisees a la Physikalisch-Technische Bundesanstalt 
(PTB) et au National Institute of Standards and Technology (NIST) ont permis d'etablir la frequence d 'horloge 
du ^''Ca parmi les frequences optiques les mieux connues (avec une exactitude relative de 1.2 • 10~^*), ces mesures 
de frequence etant en bon accord dans les deux laboratoires dans la limite de leur incertitude respective. Poten- 
tiellement, une amelioration par plusieurs ordres de grandeur de 1 'exactitude relative de I'etalon semble possible. 
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1. Introduction 

The enormous progress in laser spectroscopy, laser cooling and trapping of ions and neutral atoms, combined 
with novel efficient methods for precision optical frequency metrology has had a strong impact on the development 
of high-performance optical frequency standards and optical clocks. Such standards are widely used and serve as 
the most accurate realizations of the length unit [1], for tests of fundamental theories and for the determination 
of fundamental constants with increasing accuracy [2,3]. In contrast to today's atomic clocks operating in the 
microwave region, the optical frequencies, higher by roughly five orders of magnitude, can lead to a corresponding 
increase in frequency stability. Moreover, higher clock frequencies should lead to reduced fractional systematic 
shifts in the measured frequency, thus increasing clock accuracy. 

However, to realize the full potential of a high-performance optical frequency standard, two critical issues need 
to be resolved. First, residual line broadening and shifts caused by the thermal motion of the atoms need to be 
reduced or eliminated. Second, technical noise sources in the interrogation oscillators and the detection have to 
be suppressed. Laser-cooling techniques to achieve microkelvin temperatures have been demonstrated both for 
neutral atoms [4,5,6] in magneto-optical traps (MOTs) and single ions [7,8] trapped in the center of radio-frequency 
(Paul) traps [9]. The former systems probe a large number of atoms 10^) released from MOTs in field-free 
environments, therefore achieving exceptional signal-to-noisc ratio (SNR) but suffering from a residual Dopplor 
effect [10,11]. The single-ion approach provides a Doppler-free environment but suffers from the limitation in SNR. 
To combine the advantages of both approaches a new concept has been proposed to trap ensembles of neutral 
atoms in an array of optical dipole traps. Initial experiments based on this approach are now under way in several 
laboratories (e.g. [12]). 

Neutral atoms as well as single ions have been successfully applied in time and length metrology. Absolute 
frequency measurements with relative frequency uncertainties on the level of 10~^* have been performed for 
neutral Ca [3,13] and H atoms [14] and for Hg+, In+, Sr+, and Yb+ ions [3,15,16,17]. 

A recent advance for absolute frequency measurements in the optical domain utilizes optical frequency comb 
generators based on mode-locked femtosecond lasers [18,19]. Such a system allows one to effectively eliminate 
the uncertainty contribution of the measurement process. The capability of such systems to compare optical 
frequencies with a fractional frequency reproducibility in the range of 10~^® [20] has been shown recently. Hence 
the application of optical transitions for atomic clocks [21] becomes feasible on a level that today is only reached 
by the best atomic clocks operating in the microwave range. 

In this article, we present the results from two different realizations of an optical frequency standard based on 
ensembles of laser-cooled neutral Ca atoms released from magneto-optical traps at PTE, Braunschweig, Germany 
and NIST, Boulder, Colorado, USA. The limitation in the achievable accuracy due to the residual Doppler effect 
has been overcome on a 10"'^^ level in frequency measurements performed at PTB with a new cooling technique 
[6]. They might be minimized further to well below 10"^® by use of a new method developed at NIST based on 
an idea by Trebst et al. [10]. Frequency measurements performed in both places lead to the Ca clock transition 
frequency being known with one of the lowest uncertainties in the optical range and agree well within the respective 
uncertainty budgets of the two setups. We will summarize the results achieved for millikelvin as well as microkelvin 
atomic ensembles utilized in the frequency standards and present the prominent properties of this Ca clock 
transition that make it particularly suitable for the realization of a future optical clock. 



2. Experiment 

Both experiments use a continuously repeated measurement sequence consisting of two trapping and cooling 
stages, manipulation of the cooled atoms for investigating systematics, spectroscopy and detection. In the first 
cooling stage atoms arc captured from a thermal atomic beam and cooled to millikelvin temperatures by use 
of a standard MOT based on the broad 423 nm transition (Fig. 1) that is almost closed. At PTB, 500 mW of 
423 nm light from a frequency-doubled Ti:sapphire laser is used. A small loss channel with a small fraction of 
excited atoms on the cooling transition decaying to intermediate long lived states [22] is closed by use of a repump 
laser at 672 nm (see Fig. 1). At NIST, 40 mW of frequency-doubled light from a diode laser system is applied. 
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Figure 1. Energy level diagram of Ca with the laser transitions relevant for the Figure 2. Sketch of the MOT and spec- 

optical frequency standards at PTB and NIST. troscopy setup. 

no repump laser is used and a slowing beam is applied to increase the number of atoms in the thermal beam 
with velocities below 30 m/s that can be trapped by the 423 nm MOT. In a second cooling stage the atoms 
are further Doppler cooled in a MOT utilizing the narrow 657 nm clock transition to about 10 fiK. Here the 
excited state is optically quenched to increase the cooling and trapping forces well above the gravitational force. 
The two systems use different quenching transitions to increase the photon scattering rates (see Fig. 1). From 
the initial numbers of 8 • 10^ and 5 ■ 10® atoms at 3 mK loaded in 350 ms at PTB and 25 ms at NIST up to 
40 % are transferred to the microkelvin ensemble. The times for the second-stage cooling are 20 ms and 7 ms and 
temperatures of around 12 fiK are reached [6,23]. After the external fields are switched off releasing the atoms 
from the trap and switching on a homogenous magnetic quantization field, a spectroscopic sequence is applied 
that excites some fraction of the atoms with laser pulses resonant with the clock transition. Directly after this 
spectroscopic sequence, the normalized excitation probability of the atomic ensembles is measured. The detection 
scheme used for this purpose takes advantage of the high number of photons that can be scattered on the 423 nm 
transition [23,24,25] and gives a signal that is independent of fluctuations of atom number due to the normalization 
technique. Thus, the detection does not present a limitation to the principally achievable SNR given by quantum 
projection noise (QPN) [26] (see section 2.2). 

2.1. Interrogating the clock transition 

Prior to the interrogation of the clock transition all external fields are switched off to minimize systematic shifts 
due to electric and magnetic fields, and a homogeneous bias magnetic field of 0.1 - 0.2 mT is applied. To allow 
for the eddy currents to decay and the magnetic fields to settle, waiting times of 200 with magnetic field coils 
located in the vacuum are used at PTB, and waiting times of 1 ms are used at NIST with the coils located outside 
the vacuum chamber. To probe the clock transition at a resolution close to its natural linewidth, the use of both 
highly coherent laser radiation and atom interferometry in the time domain is necessary. 

2.1.1. Clock lasers 

The highly coherent radiation to interrogate the atoms is generated from diode laser systems whose frequencies 
are stabilized to non tunable reference Fabry-Perot resonators by means of the Pound-Drever-Hall technique 
[27]. The temperature-stabilized, vibrationally- and acoustically-isolated cavities [28], with finesses around 10^, 
provide low noise, high-Q reference resonators for the diode laser systems [24]. The laser systems are based 
on anti-refiection-coated commercial laser diodes in an extended cavity design and provide cost-effective highly 
coherent sources. With the Pound-Drever-Hall stabilization technique linewidths as low as 1 Hz (see Fig. 3) can 
be produced. The laser frequencies are widely tunable due to the application of acousto-optic modulators (AOMs) 
to provide frequency offsets between laser sources and reference resonators. Both laboratories utilize such lasers 
as highly stable master lasers. Several slave lasers are injection-locked to the master lasers in order to obtain 
the higher power (tens of milliwatts) required for the interrogation of the clock transition and the second-stage 
cooling. 

2.1.2. Atom interferometry in the time domain 

A Doppler-free method that can be described as optical Ramsey excitation [29] or as an atom interferometer 
[30] is used to overcome the limitations due to the residual Doppler broadening of about 2 MHz and 150 kHz for 

3 



1.0- 
0.8- 
0.6 
0.4 
0.2- 
O.O 





\ 

^1.5 Hz 
^FWHM 





-20 -10 10 20 
beat frequency (Hz) 




-600 -400 -200 200 400 600 
Avl (kHz) 



io-'2^ 

10-" 
10-'^ 



10 



15 



10- 



A ^ A A 



10 100 

T (S) 



1000 



Figure 3. Beat note between two diode 
lasers, locked to two independent refer- 
ence resonators (PTB). 



Figure 4. Excitation probability for the 
niicrokclvin ensemble as function of laser 
detuning for the asymmetric four-pulse 
atom interferometer with 11.6 kHz res- 
olution and pulse length set to 2.2 fis, 
obtained at NIST. 



Figure 5. Combined relative Allan stan- 
dard deviations. PTB: during frequency 
measurement of Ca-locked laser against 
H-maser (squares), cavity-locked laser 
against Ca-locked laser (triangles). NIST: 
lasers locked to Hg"*" (see [3]) and Ca, re- 
spectively (stars). 



the millikelvin and microkelvin regimes, respectively. In the time domain a sequence of two pairs of laser pulses 
from opposite directions are applied to the atomic ensemble. Each pulse provides an excitation probability of 50 % 
for atoms on resonance (Rabi-angle of 7r/2), thereby preparing each atom in a superposition of ground and excited 
states, and redirecting the atomic matter waves due to the photon recoil on two closed interferometer paths of the 
Mach-Zehnder type. For each atom passing through either of the closed interferometers the excitation probability 
at the exit depends on the accumulated phase differences between the partial waves gathered along the paths: 



^1234 = 47rr(z/L - i^Ca ± Srcc) ~ 0A(il, Ti) + 2 , - 06(^3, ^s) + 06(^4, r4). 



(1) 



The phase difference depends on the detuning of the laser frequency i/l from the clock transition vca. including 
the recoil shift Srcc = Tik'^/Aivm, with m the mass of the Ca atom, and a phase contribution seen by the atom 
due to the instantaneous laser phases (j)d at the time ti, and the atomic position r{ti) of each pulse from direction 
d = (A,B). The interference pattern, i.e., the fractional excitation as function of the laser detuning, has a cosine 
shape and is basically symmetric with its center at the clock transition (Fig. 4) when the contribution of the 
instantaneous laser phases is zero. Thus, for the realization of the frequency standard, the laser can be locked to 
the central fringe. The slight asymmetry of the envelope of the excitation spectrum shown in Fig. 4 is induced by 
atomic recoil in the absorption and stimulated emission processes and is an inherent effect. The position of the 
central fringe with respect to the frequency of the undisturbed clock transition is shifted by less than 100 mHz. 

In order to achieve almost equal excitation probabilities for all atoms in the Doppler distribution, the widths 
Tp of the squares pulses are set to 1 to 2 /xs to achieve Fourier transform widths of up to 1 MHz. By setting the 
interval T between the parallel pulses to several hundred microseconds the interferometer resolutions can be well 
below 1 kHz. 

A critical aspect of this work hinges on the observation that a non-zero offset phase shifts the fringe pattern 
with respect to vca. and leads to a systematic shift in the measured transition frequency given by 



- ( - (jiA 1 , r 1 ) + <;/iA (i2 , rz ) - 03 (ts , rs ) + (?!)B (t4 , r4 ) ) /(47rT) . 



(2) 



Quantifying the contributions to this shift is an essential part of determining the uncertainty of the frequency 
standard. At present, a direct laser light phase shift and an indirect shift due to the Doppler effect are the biggest 
contributors to the uncertainty budget (see section 3). 



2.2. Fractional frequency stability 



A key ingredient to a frequency standard is good stability. The associated short averaging times required to 
reach a given frequency uncertainty greatly accelerate measurements and help to identify systematic shifts more 
efficiently. Due to the large number of absorbers A^o, optical frequency standards with neutral atoms have the 
potential to provide a much better SNR and consequently a lower instability as compared to standards based on 
single trapped ions. The quantum-projection noise involved in the readout of the spectroscopic signal [26] sets the 
limit for the minimum fractional frequency instability achievable to l/(7r-Q-SNR), with the quality-factor Q being 
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the ratio of the transition frequency (v) to the hnewidth (~ 1/{4T)), and SNR oc Nqt /TcyrA, where t denotes 
the averaging time and Tcyci the cycle time. Using the normalized detection scheme, the detection process does not 
limit the performance of the standard. This quantum projection noise limit has already been demonstrated with a 
neutral atom microwave clock [31]. To reach this limit in an optical clock, however, extremely narrow-band laser 
sources are required. Otherwise, the Dick effect [32] will drastically limit the achievable stability. This is caused 
by aliasing of high-frequency laser frequency fluctuations introduced by sampling of the instantaneous laser offset 
from the line center in discontinuous interrogation schemes. For typical parameters of the PTB setup with the 
laser frequency noise dominated by low-frequency acoustic and seismic noise leading to a laser Hnewidth of 1 Hz 
(see Fig. 3) and the poor duty cycle between cooling and interrogation, the degradation of the stability duo to the 
optical Dick effect [33] results in a relative Allan deviation [34] a-y(r) = 2 • 10"^'* r~^^^ . At present a stability 
of C7y(100 s) = 3 • 10~^^ has been observed. With a cycle time reduced to 15 ms (but using a smaller number 
of atoms) at NIST a fractional Allan-deviation of (Ty(10 s) = 3 • 10~^® has been observed, that averages down 
for integration times of 800 seconds to 3 • 10~^* (see Fig. 5). These measurements were referenced to laser light 
from the Ifg"^ single- ion frequency standard [35] by use of a femtosecond- laser comb [20] . In both Ca systems the 
deviations from the QPN-limits are consistent with the contribution of the optical Dick effect. Even taking the low 
atom number at NIST and the Dick effect contribution into account, measurements of systematics (comparison 
of two data sets) with a relative uncertainty of 2 • 10~^® can be achieved within 200 s. 



3. Uncertainty of the calcium optical frequency standard 

In this section we present the results of frequency measurements and investigations of systematics obtained in 
both laboratories utilizing ensembles of millikelvin and microkelvin atoms. We discuss the various sources of the 
systematic shifts and their corresponding contributions to the current overall uncertainty budget. Since one of 
the primary goals of this research is to identify potential limitations for neutral atom frequency standards, we 
conclude the section with both present and projected uncertainty budgets for the Ca optical standard. 

3.1. Residual first-order Doppler shifts 

One significant contribution to the uncertainty of the Ca frequency standard is connected with residual atomic 
motion, and it is usually referred to as residual first-order Doppler effect. This is to distinguish it from the 
relativistic second-order Doppler effect due to the different inertial systems of the laboratory and the moving 
atoms. The second-order Doppler effect is negligible for laser-cooled ensembles in the millikelvin and microkelvin 
regime as opposed to thermal beams. It shifts the frequencies of millikelvin and microkelvin ensembles by less 
than 5 mHz and 20 /^Hz, respectively. 

The residual first-order Doppler effect is due to the ballistic movement of the atoms released from the MOT 
through the spatial wave fronts of the laser beams. During the time of the atom interferometry (on the order of 
1 ms) each atom "sees" a different spatial laser phase in each interaction with one of the four laser pulses that 
contribute to the overall phase shift in eq. 1. Approximating the wave fronts of the two beams at the position 
of the atomic cloud by spherical waves with wave vectors kd and radii of curvature Rd, one finds a quadratic 
dependence of the phase shift on the rms-width of the velocity distribution of the ensemble in a plane tangential 
to the curvature. A small tilt /3 of the wave vector with respect to the gravitational equipotential plane leads 
to a shift proportional to the angle and independent of the initial velocity of the atoms. Both shifts increase with 
the pulse separation T as T^. 

As has been pointed out by Trebst et al. [10], one can actually use the atoms themselves to quantify and reduce 
these shifts. By employing different types of atom interferometers, one can evaluate the laser-beam parameters 
at the location of the atomic cloud. A particularly useful interferometer uses three parallel laser pulses separated 
by a time T from a single direction. The excitation probability in this symmetric three-pulse atom interferometer 
[36,37] does not depend on the laser frequency but on the instantaneous laser phases as $123 = — 0A(ii,ri) -|- 
2</>A(t2,r2) — <j)A{t3,r3) (compare eq. 1). It has the same dependence on tilt and radius of curvature as the four- 
pulse atom interferometer. Thus it can be used to measure the contributions for either of the two beam directions 
d = A or d = B (eq. 1). In frequency measurements performed at PTB in 2001 using 3 mK atomic ensembles, these 
atom intcrforomotors were used to determine the contribution of the Doppler effect to the systematic frequency 
shifts, and a correction with a relative uncertainty below 10"^'* was applied [11]. A contribution that cannot be 
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Figure 8. Sensitivity of frequency shifts 
as function of mutual beam tilt. Launch- 
ing at 1 m/s perpendicular to the beams. 
The points are measurements for tilting 
the beams along the launch direction. 



quantified with the three-pulse atom interferometers is due to the mutual angle 9 between the counterpropagating 
beams, resulting in a non-zero sum-vector of the two wave vectors. An average velocity of the ensemble along 
this sum-vector leads to a phase shift proportional to the tilt angle. Nevertheless, by carefully overlapping the 
counterpropagating beams the uncertainty contribution with the millikelvin ensemble was reduced to well below 
10"'* [11]. 

In subsequent frequency measurements in 2003 at PTB utilizing atoms at a temperature of 12 /iK, the influence 
of the velocity-dependent wave-front curvature was drastically reduced. Fig. 6 shows this for the example of a 
particularly high wave-front curvature. Even with the very large radius of 12 m the contribution to the overall 
Doppler effect is less than 10"^"" for the microkelvin ensemble. Consequently one is left with only the contribution 
from the angle to gravity (3. In measurements at PTB, the symmetric three-pulse atom interferometer was used 
to align 13 to within 100 ^rad. This helped to reduce the total uncertainty level due to the residual Doppler effect 
to 2 ■ 10"^^ 

At NIST, this approach has been taken a step further by launching the clouds of atoms in order to enhance the 
sensitivity of the atom interferometers. Additionally, launching allows one to evaluate the various beam parameters 
independently, which enables the following straightforward optimization procedure. One first reduces the curvature 
of the wavefronts by using the three-pulse interferometer with atomic clouds launched up to 3 m/s. Due to the 
quadratic sensitivity of the phase on velocity, one can readily measure radii of curvature Rd of thousands of 
meters (see Fig. 7) with NIST's present signal-to-noise ratio. With the effects due to curved wavefronts thus 
reduced, subsequent measurements using the four-pulse interferometer can then enable reduction of the overlap 
angle 9 to less than 1 /irad (Fig. 8), corresponding to a frequency shift of 0.7 Hz at a launch velocity of 1 m/s. 
Finally, the contribution of the tilt angle can be largely suppressed with a beam-reversal technique due to the 
minimized radius of curvature and overlap angle contributions. In principle this sequence allows one to reduce the 
Doppler-shift contribution to the fractional frequency uncertainty to below 10~^^, thereby virtually removing one 
of the most challenging systematic effects for the standard. However, detailed measurements have revealed that 
the present laser-beam quality prevents performance below a level of 10~^^. While the sensitivity of the method 
and the prospects for laser beams of higher quality are certainly good, the issue of beam quality clearly remains 
one of the most critical for the future performance of the standard. 



3.2. Laser light phase shifts 



In addition to the influence of the Doppler effect on the resulting phase shifts, a direct shift of the global phase 
of the laser beams has to be taken into account. As the groups at NIST and PTB have investigated, such a phase 
shift can occur from a phase chirp produced in the switching AOMs [3,38,39]. A ringing on the edges of the RF 
pulses (and of the acoustic wave in the AOM) is imprinted on the laser pulses. This leads to a constant total 
offset phase and, hence, to a frequency shift proportional to the resolution in the four-pulse atom interferometer, 
i.e., proportional to . By measuring the frequency shift as a function of resolution and extrapolating to zero 
fringe linewidth, one can correct for this effect. Consistent results using this correction have been achieved at NIST 



6 



and PTB at respective levels of 10~ and 3- 10~ . Clearly, reduction of these uncertainties emerges as one of the 
most important issues for future frequency measurements. The determination and correction of these effects, on 
the other hand, will be limited ultimately by the remaining uncertainty of the residual first order Doppler effect, 
since both effects depend on the pulse separation T. 

3.3. Frequency shifts due to external fields 

External fields interacting with the atoms can shift the frequency of the atoms during interrogation with respect 
to unperturbed atoms. On the one hand, constant or slowly fiuctuating magnetic and electric dc fields shift the 
frequency due to magnetic effects and the quadratic Stark effect. On the other hand, thermal radiation from the 
environment and the oven as well as monochromatic radiation from the lasers can shift the frequency due to the 
ac-Stark effect. 

3.3.1. Magnetic field shifts 

The ^So {m — 0) to "^Pi (m = 0) clock transition does not have a linear Zeeman shift and has only a small 
second-order Zeeman shift. The value of the second-order Zeeman shift of the ^Pi state is known from theoretical 
calculations [40] to be (-1-63.75 ± 0.09) Hz/mT^. The diamagnetic shift of the ^So state is likely to be below 1 % 
of this value. Measurements on atomic beams and laser-cooled atoms [41,42,43] show the second-order Zeeman 
shift of the clock transition to be (-1-64 ± 1) Hz/mT^. Typically a quantization field of 0.1 to 0.2 mT separates 
the Zeeman components (Am — ±1) well from the Am = component. This field can be measured by use of the 
Zeeman components that are separated from the Am = component with ±21.01 MHz/mT [41]. When measuring 
the field with a relative uncertainty of 0.5 % and using the theoretical coefficient with an uncertainty of 1 %, the 
quadratic Zeeman shift at 0.2 mT can be corrected with an uncertainty of 36 mHz. 

3.3.2. DC electric fields 

The coefficient for the quadratic Stark effect is well known from measurements by Zeiske [44] of the differences 
Aa of the static polarizabilitics of the clock states. For electric stray fields stemming from arbitrary electric 
charges on the surfaces of the vacuum chambers and magnetic field coils (distance to trap > 3 cm) an average 
value Aa = (3.07 ± 0.06) //Hz/(V/m)^ can be used to correct for shifts. If we assume that electric potentials are 
less than 10 V on the grounded surfaces {E < 200 V/m), the resulting shift lies between mHz and -70 mHz. 
If we assume equal probability of a shift within this interval, the average shift is (—35 ± 20) mHz. Higher-order 
contributions shift the transition by not more than a few microhertz for dc field levels up to 1 kV/m. 

3.3.3. Blackbody radiation 

The infiuence of the radiation due to the ambient temperature in the lab (between 288 K and 298 K) can be 
estimated based on a quasi-static approach [44]. With the atoms placed in a blackbody- radiation environment 

of temperature •& the systematic shift is (—1.31 ± 0.03) • 10^^" ^| ■ i?*. The infiuence of the radiation of the Ca 
oven at 900 K can be obtained firom an exact calculation based on an approach by Farley and Wing [45], which 
includes all relevant transitions. The calculated shift is about 45 % larger than the one found with the quasi-static 
approach. However, it is difficult to estimate the level of radiation from the oven reaching the atomic cloud. 
Therefore measurements of the shifts at various oven temperatures were performed at PTB during the frequency 
measurements in 2001 and 2003, resulting in an estimate of the shift of (—2.8 ± 3.9) Hz. This is at present the 
largest contribution to the total uncertainty. Proper shielding of the trap from the oven, possible with the use of a 
Zeeman slower and an optically deflected atomic beam, will allow ont; to eliminate the oven contribution. This will 
reduce the uncertainty in the correction of the blackbody shift to the contribution of the room-temperature. For 
a temperature of i? = 293 K known to ±1 K the uncertainty in the correction could then be reduced to 24 mHz. 

3.3.4. AC Stark effect 

In the PTB setup the infiuence of residual laser light from the cooling and quenching lasers is eliminated by the 
use of mechanical shutters. At NIST the light is switched off using multiple deflections in AOMs. Measurements 
with increased residual light levels, when AOMs were left on, showed that even when the frequency offsets with 
the AOMs switched off are neglected, the shifts are loss than 1 mHz for the blue cooling laser and 85 mHz for 
the quench laser. Presently, the red light in the cooling and probing beams is less well isolated and can lead to 
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Figure 9. Measurement of the change of 
the Ca clock frequency with a chang- 
ing flux of thermal atoms. Accessed 
via NIST's femtosecond-Iascr comb, the 
mercury single-ion frequency standard 
was used as a reference. The difference 
found for doubling the atomic flux is 
(-0.9 ± 2.2) Hz. 



Figure 10. Measurement of the density-dependent frequency shift (a) using cold 
atoms {T ^ 3 mK), (b) using ultracold atoms (T ^ 20 /xK). The lines represent a 
linear regression to the data 



frequency shifts of up to 0.5 Hz. Cascading multiple AOMs in a suitable way will allow one to reduce this shift 
to the millihertz level as well. Hence, even at a relative uncertainty level of 10"^^, the use of mechanical shutters 
can be avoided, thereby avoiding compromises regarding cycle times and stability. 



3.4. Collisions with thermal atoms 



In the PTE setup the cold ensemble is shielded from collisions with thermal atoms by a beam block placed in the 
path of the thermal beam. Without the aperture, trap losses due to collisions with thermal atoms would strongly 
reduce the number of atoms in the trap. In the similar NIST setup no aperture is used, but the oven is run at a 
lower temperature. The influence of collisions with thermal atoms on the frequency shift was measured referencing 
the Ca frequency to NIST's optical Hg""" frequency standard and to NIST's H masers using the femtosecond-laser 
comb. Fig. 9 shows a measurement against Hg""" when the flux of atoms was changed by a factor of 2. No frequency 
change was observed that would be significant on a 1 Hz level. In the future it can be completely avoided by using 
a deflected slow atomic beam for loading the MOT. 



3.5. Influence of cold and ultracold atomic collisions 



Interactions between the atoms in the laser-cooled ensemble, i.e., cold collisions, lead to a frequency shift 
of the clock transition. The magnitude of the shift depends on the states of the two colliding atoms and their 
separation, and results in a mean frequency shift that depends on the temperature and the density of the atomic 
ensemble. For the low densities ( about 10^'' m~^) in PTE's ensemble, the frequency shift is assumed to be a 
linear function of the collision rate, and hence of the density p. In order to determine the density dependence 
of the frequency shift, the frequency differences were measured as a function of average density differences [11]. 
The density was varied by switching the repump laser at 672 nm (Fig. 1) on or off during the loading time of the 
millikelvin ensemble. This changed the number of trapped atoms while all other parameters such as the magnetic 
fleld, loading time, frequency of the trapping light and laser-beam proflle were kept constant. The radius and 
temperature of the trapped cloud were checked and did not change signiflcantly when the average density Ap 
was varied by a factor of five. The measurement of the coUisional shift Ai/fi/ depending on the density p (Fig. 
3.4 a) gave Av/u = (3 ±4.4) • 10"'^° m^ x p [11]. The corresponding frequency shift for the average density of the 
expanding ensemble when measuring at a resolution of 1.2 kHz (T = 215 ps) would be (1.2±1.8) Hz. Hence, within 
the uncertainty of this measurement, no significant shift has been observed. The same was true when microkelvin 
atoms were used (Fig. 3.4b) where an even lower limit of Av jv — (—0.4 ± 1.2) ■ 10"''^ m^ x p was obtained. At 
the usual densities at the PTE setup this would lead to a shift of (20 ± 60) mHz, while in the NIST setup the 
densities and therefore also the shifts are lower by at least one order of magnitude. Within its uncertainty this 
coefficient is four orders of magnitude smaller than that seen in Cs-fountain clocks [46], which have nonetheless 
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Uncertainty (Hz) 


Projected Uncertainty (mHz) 




NIST 


PTB 


PTB 






Oct. /Nov. 2000 


Ocl. 2001 


Oct. 2003 




Effect 


(T = 3 mK) 


('/' = 3 mK) ('J' = 12 pK) 


(T = 12 pK) 


Ist-order Doppler effect 


12.1 


2.6 


1.0 


4 


2nd-order Doppler effect 


0.005 


0.005 


2 • 10"^ 


0.02 


other phase contributions 


4 


4 


1.6 


5 


asymmetry of line shape 


0.05 


0.05 


0.05 


5 


2nd-order Zeeman effect 


1 


0.1 


0.1 


36 


AC-Stark effect 


2 


0.1 


0.1 


1 


quadratic Stark effect 


0.06 


0.02 


0.02 


20 


blackbody radiation from 










oven 


1 


4.3 


3.9 


- 


cliambcr 


0.07 


0.07 


0.07 


24 


collisions of cold atoms 


10 


1.8 


0.06 


10 


electronic stabilization & laser drift 


10 


3.2 


0.1 


1 


sum in quadrature 


19.2 


7.4 


4.3 




statistical uncertainty of frequency mcas. 


2.5 


3 


3 


0.05 


Cs clock (about 1 ■ 10"^^) 


0.8 


0.5 


0.5 




total uncertainty 


19.4 


8 


5.3 




total relative uncertainty bv jv 


4.3 • 10"" 


1.8- 10-^* 


1.2- 10-" 





Table 1 

Uncertainty budget for frequency measurements utilizing the Ca frequency standards at PTB and NIST with millikelvin atoms 
(col. 1, 2) [25,38], and with microkelvin atoms (col. 3). Column 5 summarizes the projected uncertainty budget for an optimized 
Ca frequency standard (see text). 

achieved an uncertainty level of 10~^^. Thus one sees the potential of this optical frequency standard to achieve 
very low uncertainty levels while maintaining a very high stability with a high number of atoms lO'^). 

3.6. Frequency measurements and uncertainty budget 

Table 1 summarizes the results for the present uncertainty budget for frequency mcEisurements with millikelvin 
atoms in 2001 at PTB [47,13] and 2000 at NIST [38,48,3], as well as a recent frequency measurement with 
microkelvin atoms in 2003 at PTB. While the measurements at NIST in 2000 were still limited by the unknown 
effect of cold collisions and the exact contributions by the Doppler effect, the PTB measurements of 2001 were 
limited by the blackbody shift (4.3 Hz) duo to the oven being close to the atomic ensemble. This contribution 
still dominates the uncertainty budget in PTB's frequency measurement of 2003, while the Doppler effect and 
subsequently other phase contributions were further reduced by use of microkelvin atomic ensembles. A setup with 
the oven removed from the vicinity of the trapped ensemble is currently being assembled. Asymmetries in the 
excitation spectra that can lead to systematic frequency shifts when stabilizing the laser have been estimated as 
well as measured [42] indicating a shift of less than 50 mHz with a 3f-stabilization scheme. The use of femtosecond- 
laser frequency combs for the comparison against state-of-the-art Cs-fountain atomic clocks (and against the Hg"*" 
single-ion frequency standard at NIST) allowed both groups to avoid being limited by either the primary time and 
frequency standards or the statistical uncertainty of the comparison. The results are summarized in Fig. 11 together 
with the first frequency measurements of the Ca clock transition performed at the PTB, using phase- coherent 
frequency multiplication chains [49,42]. The average over all measurements is i'ca = (455 986 240 494150.2 ± 
7.7) Hz. The latest measurement with microkelvin atoms performed at PTB in October 2003 resulted in vca = 
(455 986 240 494 144 ± 5.3) Hz. The results of all measurements agree well with the average value within their 
respective uncertainties (sec fig. 11). This makes the Ca optical frequency standard with neutral atoms the best 
such standard evaluated in more than one laboratory. 
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Figure 11. Results of frequency measurements of the 657 nm Ca intercombination line performed at PTB (solid circles) and NIST 
(open diamonds) with the weighted average (solid line) of all measurements and the statistical uncertainty (dashed lines). 



3.7. Projected uncertainty budget 



The projected uncertainty budget in the last column of Table 1 is based on the present knowledge on the 
system as presented in the above sections. Possibly, the dominating uncertainty contribution will result from the 
correction needed for the room-temperature blackbody shift. The limit is due to the uncertainty in the correction. 
Modelling of the atomic signal (see [25]) based on a description by Borde et al. [29] to correct for shifts due to 
the asymmetries in the line shape (see Fig. 4) will allow one to quantify the shift more accurately. A correction 
of the shift to within 2 % would reduce the uncertainty to 1 mHz. 

Electronic stabilization to line center at a level better than 10~® has already been demonstrated in the NIST-7 
thermal atomic Cs-clock [50]. An offset from the line center due to drift of up to 10 Hz/s of the free-running laser 
can be reduced to below one millihertz when two integrator stages and measurement cycles on the 10 ms range 
are used (see [42] for a mathematical model describing such a stabilization scheme). 

With 10^ atoms confined to a trap of radius 2 mm the present knowledge on density shifts would be sufficient 
to make this uncertainty contribution negligible. Furthermore, local higher-order distortions of wave fronts on a 
scale of 1 mm that cannot be clearly resolved by launching the ensemble would be averaged away over the size of 
the trap. A 12 ensemble with an average velocity of less than 1 mm/s would enable one to achieve a negligible 
residual first-order Doppler effect. Nevertheless a wave-front flatness with Rd > 3000 m would need to be achieved 
simultaneously and would certainly present a technical challenge. 

The numbers given in the table are best-case estimates based on our present knowledge of what could be 
achieved. Whether all these conditions can be simultaneously achieved, as would be required for a total combined 
uncertainty below 10"^^, is subject to further investigation. 



4. Conclusion 



We have described the current status of the optical Ca frequency standard with laser-cooled neutral atoms as 
realized in two different laboratories for the purpose of developing a possible future optical atomic clock. Fre- 
quency measurements performed make the frequency of the clock transition of ''"Ca one of the best known optical 
frequencies (uncertainty 1.2 • lO"^**), and agree in both laboratories to well within their respective uncertainties. 
Furthermore, extensive investigations of systematics carried out in both laboratories do indicate that the presently 
visible roadblocks are of a purely technical nature and that measurement capabilities to identify these shifts al- 
ready exist. With the present knowledge of the properties of this frequency standard an attainable uncertainty 
level of below 10~^^ can be anticipated. 

Taking advantage of the unprecedented potential stability of 10^^^ in 1 s that is feasible with this standard, the 
systematic uncertainty level could be reached in a few seconds of measurement time. This would however require 
a free running spectroscopy laser of low frequency noise during the 100 ms needed for a full stabilization cycle. In 
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case of white frequency noise this would correspond to a laser linewidth of less than 100 mHz. 

Such a standard is a promising candidate not only for an application as a future atomic clock, but as well as a 
versatile instrument for precision measurements including the search for variation of the fundamental constants 
and the test of fundamental theories (e.g. see [51] and references therein). 
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